Dynamic concentration of motors in microtubule arrays 
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We present experimental and theoretical studies of the dynamics of molecular motors in microtubule 
arrays and asters. By solving a convection-diffusion equation we find that the density profile of 
motors in a two-dimensional aster is characterized by continuously varying exponents. Simulations 
are used to verify the assumptions of the continuum model. We observe the concentration profiles 
of kinesin moving in quasi two-dimensional artificial asters by fluorescent microscopy and compare 
with our theoretical results. 
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The cytoskeleton is a network of polymers essential for 
the dynamic organization of many eukaryotic cells. Its 
function depends not only on protein fibers, but also on 
many accessory components [jjj . Among these are motor 
proteins that reversibly bind to and walk along the sur- 
face of cytoskeletal polymers, consuming ATP as a source 
of energy. A natural consequence of this directed move- 
ment in organized fiber arrays is an non-uniform spatial 
distribution of the motors. In vivo microtubules are of- 
ten observed in radial arrays, or asters, where all mi- 
crotubule "minus" ends are at the center, and the "plus" 
ends are radiating outward. Aster of opposite polarity, 
in which kinesin moves inward can also be formed in- 
vitro Q. 

In this letter we present quantitative experiments on 
a quasi two dimensional aster, fig. (|l|), formed from 
long microtubules polymerized between closely separated 
cover slips. We analyze the results of these experiments 
using a 2D transport diffusion equation. The theory con- 
tains a number of approximations, due to angular averag- 
ing and projection from 3 to 2 dimensions. We verify that 
no major quantitative error is introduced by perform- 
ing simulations in the full confined geometry. We also 
consider, theoretically, the case of one and three dimen- 
sions: The aster in one dimension corresponds to a tube 
in which microtubules are all oriented in the same direc- 
tion. This is the case for example in the axons of nerve 
cells. As distinct motors move toward the "plus" or the 
"minus" end, we consider both cases of inward/outward 
directed motion. 

Consider TV immobile straight microtubules radially 
arranged in the available volume. Molecular motors are 
present which can exist in two different states, either at- 
tached to a filament, or detached. Unattached motors 
diffuse freely, with a diffusion constant D. Attached mo- 
tors move on their filament (radially in the aster geome- 
try) at a velocity v. Positive values of v corresponds to 
outward movement. Transitions between the two states 
are stochastic: The motor spontaneously detaches from 
the microtubule at an unbinding rate p oS (s _1 ). Far from 




FIG. 1. In the presence of a microtubule array motors 
can move by free diffusion in solution or by directed motion 
on microtubules. Movement of motors in an aster can lead 
to accumulation, if the motor moves inward (top right), or 
depletion, if the motor moves outward (bottom right). Accu- 
mulation also occurs in oriented parallel microtubule arrays 
(bottom left). 

saturation, the number of binding events per second is 
proportional to the local concentration of free motors, 
and to the number of available binding sites on the micro- 
tubules. If the concentration of free motors is expressed 
in molecules per cubic micrometers, and the available 
"quantity" of microtubules in micrometers, the constant 
of proportionality p on has the dimension of a diffusion 
constant ||,[|. 

Let b(r) and f(r) be the concentrations of bound and 
free motors, respectively, at distance r from the cen- 
ter, averaged over all angles. Bound motors move radi- 
ally at speed v, and create a convective flux J5 = vb. 
Unbound motors diffuse freely, creating a radial flux 
Jf — —Ddf/dr. If S is the surface area at distance r, 
there are p°^Sb dr release events in the volume between 
the radii [r, r + dr], and p on Nf dr attachments per sec- 
ond. We therefore obtain the coupled kinetic equations: 
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To find the steady state, we use the fact that the net flux 
Jb + Jf is zero, implying b = Dv~ 1 (df /dr). Substituting 
this result into equation (Q), and denoting /' = df /dr 
and f" = d 2 f/dr 2 , we find in d dimensions 
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The physical parameters in the problem have been re- 
duced to 

a = v/p oS 

/3i = S 1 D/p on N, f3 2 =2irhD/Np on , fa = 4nD/p on N 
7 = D/v. (3) 

where Si is the area of the tube in one dimension and 
h is the sample thickness in a quasi two-dimensional ge- 
ometry, a is the average distance that a motor moves 
on a microtubule, before detaching. (3d characterizes the 
geometry of the aster, and 7 determines the relative con- 
centration of bound to free motors. Note that for inward 
movement, a and 7 have negative values, while for out- 
ward movement all parameters are positive. 

For the motor protein kinesin, experimental data pro- 
vide values for the parameters in the model. The walk- 
ing speed of kinesin without load is v = 0.8 /ims -1 
The unbinding rate p oS is obtained from the average dis- 
tance that kinesin moves before detaching. Measured 
average run length, a = v/p oS , are for kinesin in the 
range 0.4 — 1.5 /im ; we use p oE = 1 s _1 . Direct 

chemical measurements of p off j| agree with this value. 
The binding rate p on of the kinesin construct used in our 
experiment has not been directly measured. To estimate 
it, we assume that interaction between microtubules and 
motors is diffusion limited [|l0|| : Measurements with ki- 
nesin' s soluble dimeric motor domain Q provide a value 
of p on k in = 7.3 /im 2 s _1 , and its diffusion constant is 
50 /im 2 s _1 [pT| , p!2| . For single kinesin adsorbed on beads 
the equilibrium constant for the binding convection 
(equal to p on /p off ), provides p on bead = 0.25 /im 2 s _1 , 
and a diffusion constant of 2 /im 2 s _1 |l3). The ratios 
p on /D are 0.14 and 0.12 respectively, and we use the 
averaged value p on /D = 0.13. Finally, based on its 
molecular weight, we estimated a diffusion constant of 
D = 20 /iii^s -1 [@, which yields p on = 2.6 /rni 2 s _1 . 

In a tubular, quasi one dimensional geometry, all mi- 
crotubules are oriented in the same direction, (Fig [I]), 
lower left. We find that the steady state profile for the 
motor concentration is exponential f(r) ~ e r l a where the 
distance a is a root of the equation a 2 + a/a — = 0. 
To relate this result to the situation of motors in a cell, 



we now consider a tube filled with oriented microtubules 
connected on one side to a large body. Motors have a 
concentration which varies exponentially with the dis- 
tance from the cell body. The concentration at the end 
of the tube is e L l a times smaller (or greater depending 
on the motor sense) than in the cell body, where L is 
the length of the tube. We estimated a for a cellular ex- 
tension of a diameter of 2 /im, and for parameters of the 
motor kinesin: If the extension contains 20 microtubules, 
then a = 2 /im; for the same tube containing a single mi- 
crotubule, a = 30 /im. Therefore an unregulated kinesin 
(which can always bind and move) would be concentrated 
even in short extensions of a cell containing outward po- 
larized microtubules. This is indeed observed in vivo for 
kinesin heavy chain if it is over-expressed in the absence 
of the regulatory light chain |15|. 

Our experiments were performed in a quasi 2D ge- 
ometry. We thus give particular weight to the theo- 
retical analysis of this case which also presents some 
interesting theoretical features: Consider firstly a non- 
motile binding protein (v = 0). At equilibrium, the un- 
bound molecules are evenly distributed throughout the 
volume, f(r) = const , while the concentration of bound 
molecules is proportional to the local concentration of mi- 
crotubule so that b(r) ~ 1/r . Thus binding of motors 
induces their accumulation in the center, where micro- 
tubules are more concentrated. 

For general speed v the solutions of the equations 
in 2D can be expressed in terms of Whitaker func- 
tions, however, simple analysis (performed by substitut- 
ing f(r) ~ r~ x in eq. (2)) shows that the solution in 
the quasi-two dimensional case is well approximated by 
power laws beyond the radius o? j f(r) ~ r a ^ 2 , and 
b(r) ~ jf(r)/r . Thus the concentration profile of mo- 
tors is characterized by an exponent which is a continu- 
ous function of the physical parameters. From the above 
expressions we find a = ±0.8/ira, 1/ira < P2 < 10/ira, 
and 7 = ±60/im. The theory shows that in large asters 
(for kinesin, and a sample thickness of 9 /im, asters with 
more than ~ 600 microtubules), most of the motors are 
trapped in the center, and very few motors are left else- 
where in the sample, effectively causing a dynamic "lo- 
calization" . The depletion from the aster center of a ki- 
nesin motor moving outwards is comparatively weaker. 
For small asters, motors concentration can be higher in 
the center, merely as a consequence of the binding of mo- 
tors to microtubules. Total depletion is achieved only for 
large asters (of 1000 microtubules), for which outward 
transport overcomes the pure binding effect. 

In three dimensions, the perturbation in the concen- 
tration due to the presence of the aster is significant only 
within a distance a/vT^ from the center. For large radii, 
f(r) ~ (1 — r*o/r). We did not study this situation further 
due to the absence of experimental results. 

To derive the convection-diffusion equations, we aver- 
aged over the directions transverse to the microtubules 
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FIG. 2. Simulated bound (triangles), free (circles) 
and total=bound+free (squares) motor concentrations and 
2D-theory (lines), for kinesin moving inward (left) or out- 
ward (right) in asters having either 300 (top) or 600 (bottom) 
microtubules. As expected some deviations are seen near the 
center of the aster due to finite sample thickness implemented 
in the simulation. 



(both angularly, and over the thickness of the sample). 
This approximation breaks down experimentally in a thin 
sample at the center of a aster where the geometry is 
three dimensional, or at large radial distances where the 
microtubules are too far apart for angular averaging. We 
thus performed simulations in a true confined geometry 
in three dimensions to check that no substantial errors 
are introduced in the two dimensional theoretical descrip- 
tion. In our simulations the aster is formed by micro- 
tubules of length L = 50 //m, with their "plus" end in the 
center of a cylindrical box of radius L, and of thickness 
9 /im. Each motor is characterized by its state (bound 
or free), and a vector (position). To compute the bind- 
ing of motors, a rate p + and an interaction range e were 
introduced, so that at each step, a free motor has a prob- 
ability p + dt to bind to any filament located at a distance 
(by projection) closer than e, this effectively corresponds 



to p° 



At each time step dt = 4.10 s, bound 



motors may detach with a probability p oS dt, and other- 
wise move radially by a distance v dt; free motors make 
random steps, with < dx 2 >= 3D dt. The motor param- 
eters were taken to mimic kinesin (see above), with the 
additional value e = 50 nm [[|, and p + = 312 s _1 (which 
yield the correct value for p on ). Other choices of e and 
p + conserving p on gave similar curves. The agreement 
between simulation and theory, Fig. confirms that 
our analytical approximations are faithful to the model. 

We now turn to our experimental observations on con- 
fined quasi two dimensional samples: Using fluorescent 
microscopy we measured the kinesin distribution in asters 
with the two possible polarities. These asters have either 




Motor and microtubule distribution in experi- 
mental asters with different number of microtubules: Fluo- 
rescence images (90x70 /im) obtained for the motors (top), 
for the microtubules (middle). The cross mark the computed 
measured center. The intensities of the images are here scaled 
by different factors. 



the plus or the minus ends of the microtubules in the 
center [[l4|]2|, and take about 30 minutes to form. All 
data presented here are extracted from two identically 
prepared samples, in which many asters of various size 
formed. We measured 115 regular asters with an auto- 
matic epi-fiuorescence microscopic setup (Zeiss axioplan 
2 with Olympus 100X oil-immersion objective). We de- 
tected the motors (labeled with the fiuorophore fluores- 
cein) and the microtubules (labeled with rhodamine) in- 
dependently. Digital pictures were taken with a 12-bit 
CCD camera (Hamamatsu C4742-95, 1280x1024 pixels). 
The camera is linear, and unsaturated pixel values reflect 
the relative quantity of protein in the imaged region. The 
sample thickness was 9 /im. 

The center of the aster and the profiles of fluorescence 
intensity are calculated from the image. A common back- 
ground pixel value was subtracted from the motor pro- 
files, which are then normalized. An exponent is obtained 
by fitting the profile in the range 1.5-20 /im (the data 
below 1.5 /im is noisy). To measure the number of mi- 
crotubules in the aster, we fit the profile of microtubule 
fluorescence to the function (M/r + B), where r is the 
distance from the center. B is a background, and M is 
proportional to the number of microtubules. Calibration 
was done by manually counting the microtubules in five 
asters. The 1/r profile corresponds to an homogeneous 
aster of long microtubules. Experimentally the asters are 
not perfect (some, like Fig ||, left are not well focused). 
When the fit of the microtubule-profile to 1/r is poor, 
we have no reason to expect the theory to apply. These 
asters are plotted with a different symbol in Fig. (||). 

We measured the distribution of kinesin in asters con- 
taining different numbers of microtubules. Three typical 
examples are given in Fig. ( ||). Motor profiles of indi- 
vidual asters are rather well fitted by a power law. They 
are almost linear on a log- log plot, and steeper for bigger 
asters Fig. ( ||), inset. Plotting the exponent of the mo- 
tor profiles as a function of the number of microtubules 
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FIG. 4. The effective exponent of the kinesin concentration 
profile becomes more negative with the increasing numbers of 
microtubules in the aster (continuous line: 2D-theory). Filled 
symbols correspond to regular asters for which the radial mi- 
crotubule density falls of as 1/r. Open symbols correspond to 
irregular asters where the density is inconsistent with a law in 
1/r. Arrow heads points to asters shown in previous figures. 
Inset (log- log): Experimental (circles) motor concentration 
profiles, as extracted from the pictures shown in Fig. ^, and 
2D-theoretical curves computed for the measured number of 
microtubules (27, 148, 231) in the aster. 



extracted from the microtubule profiles allows us to com- 
pare directly experiments and theory (see Fig. |]). The 
data points, each representing one aster, are scattered 
around the theoretical curve, reflecting the heterogeneity 
of the asters. However, the major trend in the exponent 
is correctly predicted by the theory, so that denser asters 
are characterized by a larger localization exponent for 
the motors. We also imaged kinesin moving outwards in 
asters of normal polarity but the signal was too dim 
to extract a reliable profile (the predicted kinesin profile 
in this situation is rather flat). 

Accumulation or expulsion of molecular motors in 
asters may have important functional implications in bi- 
ology. For instance within a spindle made of two inter- 
acting asters of microtubules, minus-ended motors could 
concentrate at the poles while plus-ended motors would 
be excluded from the same regions. This could contribute 
to the mechanism of spindle assembly, and/or to its me- 
chanical stability. Interestingly we find that kinesin accu- 
mulates in asters of 300-1000 microtubules, which is com- 
parable to the number of microtubules present in spindle 
asters of most animal cells. However, the geometry and 
motors of the spindle are not the one studied here. 

We did not consider the regulation of motor activity: 
In our study motors can always bind and move on fila- 
ments. Cells use a variety of processes to counterbalance 
the impact of motor transport on their localization. For 
example, the folding of kinesin into a non-motile confor- 



mation, in the absence of a cargo [pT.17,1^] dampens the 
transport-induced localization. The recombinant kinesin 
fraction used in our experiment lacks this capacity. Even 
with this partial inhibition, the movement of loaded ki- 
nesin brings them to places from which they have to be 
recycled. Additional regulation mechanisms include local 
synthesis and degradation of the motors, involvement of 
motors of different directionality transporting each other, 
etc. On the other hand, we can also imagine situations in 
which the unregulated localization of a motors resulting 
from their movement can have interesting consequences. 

In summary, motor movements on microtubules can 
effectively cause their "compartmentalization" . The the- 
ory provides a full understanding of the influence of all 
motor kinetic parameters, and of the geometric proper- 
ties of the microtubule array. 

We thank A. Ajdari, S. Blandin, A. Desai, E. Karsenti and 
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